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A stirring device consisting of a periodic motion of rods induces a mapping of the 
£^ fluid domain to itself, which can be regarded as a homeomorphism of a punctured 

surface. Having the rods undergo a topologically-complex motion guarantees at least 
i_i a minimum amount of stretching of material lines, which is important for chaotic 

^ mixing. We use topological considerations to describe the nature of the injection 

of unmixed material into a central mixing region, which takes place at injection 
cusps. A topological index formula allow us to predict the possible types of unstable 
l/-} foliations that can arise for a fixed number of rods. 
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By stirring a fluid, mixing is greatly enhanced. By this we mean that if our 
goal is to homogenise the concentration of a substance, such as milk in a teacup, 
then a spoon is an effective way to mix. But in many industrial applications, 
such as food and polymer processing, the fluid is very viscous, so that stirring 
is difficult and costly. Hence, insight into the types of stirring that lead to good 
mixing is valuable. We explain how the mixing pattern — the characteristic 
shape traced out by a blob of dye after a few stirring periods — is tightly 
connected with topological properties of the stirring motion. In particular, we 
can enumerate the allowable number of pathways where material gets injected 
into the mixing region, as a function of the number of stirring rods. 
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FIG. 1: The figure-eight stirring protocol. The inset shows the sequence of rod motions. The central 
black circle is the stirring rod, and the other two circles locate the position of regular islands, which 
serve as 'ghost rods.' (Experiments by E. Gouillart and O. Dauchot, CEA Saclay [1].) 

I. STIRRING WITH RODS 

A rod stirring device, where a number of rods are moved around in a fluid, is the most 
natural and intuitive method of stirring. The number of rods, their shape, and the nature of 
their motion constitute a stirring protocol. For example, Fig. [I] shows the result of stirring 
with the figure-eight protocol, whereby a single rod in a closed vessel traces a lemniscate 
shape. The Reynolds number is very small, so that the fluid (sugar syrup) is in the Stokes 
regime, where inertial forces are negligible and pressure and viscous forces are in balance. 
Because of the shape of the rod and container, three-dimensional effects are negligible. The 
fluid is the pale background, and a blob of black ink has been stretched by a few periods of the 
rod motion. The evident filamentation of the blob is characteristic of chaotic advection [21 E] , 
which greatly enhances mixing effectiveness in viscous flows. 

We aim to understand the features of rod stirring protocols such as those depicted in 
Figs. [TJ |2j and [3] from topological considerations. The new topological features that we 
discuss can be divided in two broad categories: (i) The injection cusps and their role; (ii) 
The identification of higher-pronged singularities directly in flow simulations. Though both 
of these concepts are familiar from the topological study of surface homeomorphisms (H 
[5], we interpret them here in light of practical stirring protocols. Our study is a natural 
continuation of the original investigation by Boyland et al. [6] and the subsequent work of 
Gouillart et al. Our goal is to refine the previous approaches by looking for more detailed 
features of the Thurston-Nielsen classification in real fluid flows. 

As an illustration, we examine the topological features visible in Fig. [TJ The two unmixed 
regions in the centre of the loops of the figure-eight play the role of two extra rods, called 
ghost rods [7], so that we can regard this protocol as effectively involving three rods. Three 
rods is the minimum needed to guarantee exponential stretching of material lines [SI IE], and 



(a) 



(b) 



FIG. 2: Numerical simulations of flow in a channel with a figure-eight rod stirring protocol, (a) 
Injection cusp against flow; (b) Injection cusp facing flow. The flow direction is from the bottom 
to the top of the figures. 



for the figure-eight protocol the length of material lines grows by at least a factor of (l + v^) 2 
at each period. This is the first feature that can be understood from the topology of the 
rod motion: it places a lower bound on the topological entropy, which is closely related 
to the rate of stretching of material lines in two dimensions. This aspect has been well 

studied jsiciisiisiiiniiniiiaiisiiiaiiaiisiiiz!. 

Less well studied is another crucial feature obvious in Fig. [lj the unmixed fluid (white) 
is injected into the mixing region (kidney-shaped darker region) from the top part of the 
region, where a cusp is clearly visible. We say that the figure-eight protocol has one injection 
cusp. In fact, for three rods the situation in Fig. [I] is typical; for instance, a kidney-shaped 
mixing region is also evident in the efficient stirring protocol of Boyland et al. j6]. The 
nature of the injection of unmixed material into the central region has a profound impact 
on mixing rates, as was shown in recent experiments (lj. In these experiments the rate of 
injection of unmixed material into the central region of a stirring device dramatically limited 
the efficiency of mixing. This injection took place along injection cusps as described here, 
hence the number of injection cusps and their positions are clearly important for mixing. 

The importance of injection cusps is even more apparent when dealing with open flows. 
Open flows, as opposed to flows in closed vessels, involve fluid that enters and subsequently 
exits a mixing region. This situtation is very common in industrial settings, since it allows 
continuous operation without the need to empty the vessel. Typical fluid particles only 
remain in the mixing region for a short time. Figure [2] shows the figure-eight stirring pro- 
tocols in an open channel. The only difference between cases (a) and (b) is the direction 
of rotation of the rods. In a closed vessel, reversing the direction of rotation merely moves 
the injection cusp from top to bottom, but in an open flow it moves the injection cusp 



either facing the flow or against it. In Fig. 2(a), the injection cusp is downstream, so the 



passive scalar enters the mixing region from behind. This is clear from the dye filament at 
the bottom, which skirts the mixing region from the right. In contrast, Fig. 2(b) shows the 
opposite case, where the injection cusp is upstream. In that case the dye is drawn directly 
into the mixing region. Clearly these two cases have radically different mixing properties, as 
evidenced by the different dye patterns downstream. Thus, being able to predict the number 
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and nature of injection cusps from the rod motion is of crucial importance. The qualitative 
details of Fig. [2] are independent of the manner in which dye is injected. The extra injection 
cusps visible downstream in Fig. 2(a) are an artifact of the open-flow configuration: they 
are images of the single injection cusp being advected downstream by the mean flow at each 
period. The results presented in the rest of the paper apply rigorously only to closed flows, 
but it is evident from Fig. [2] that many qualitative features carry over from closed to open 
flows. 

In this paper we will see that the number and position of injection cusps depends on the 
topology of the rod motion. The number of injection cusps will depend on the number of 
rods (or ghost rods), but for a fixed number of rods only a few configurations are possible. 
This is because the number of injection cusps is constrained by a topological index formula. 
Jana et al. [TS] studied the impact of the topology of streamlines on chaotic advection. They 
used the Euler-Poincare-Hopf formula (see Section III) to determine the allowable fixed- 
point structure of steady velocity fields. When the flow is time-dependent, instantaneous 
fixed points of the velocity field mean little for chaotic advection (except when the time 
dependence is weak). Our study addresses arbitrary time-periodic flows by studying the 
mapping of fluid elements (the Lagrangian map) directly, from a topological perpective. 

The paper is organised as follows. In Section [IT] we present some necessary mathemat- 
ical background, in particular the idea of pseudo-Anosov stirring protocols. We identify 
mathematical objects such as the unstable foliation in a specific fluid-dynamical example. 
The unstable foliation associated with a pseudo-Anosov protocol is the central object of our 
study. In Section III we use a topological index formula to enumerate the possible unstable 
foliations for a given number of stirring rods, and Section IV is devoted to a stirring device 
that exhibits a hyperbolic injection cusp, that is a cusp associated with the unstable manifold 
of a hyperbolic orbit. Such a cusp is less obviously identifiable from flow visualization, but 
can still be inferred by examining the topological properties of the stirring protocol. Finally, 
we summarise and discuss our work in Section |V] 



II. PSEUDO-ANOSOV STIRRING PROTOCOLS 

In this section we will translate the physical system — a stirring protocol — into objects 
suitable for mathematical study. As mentioned in Section [I], our focus will be on very viscous 
flows, where stirring and mixing is challenging, and we consider the flow to be essentially 
two-dimensional, as is typical of shallow or stratified flows. In these circumstances, a periodic 
motion of rods leads to a periodic velocity field. (Note that the rods may end up permuted 
amongst themselves at the end of each period.) This velocity field will induce a motion of 
fluid elements from their position at the beginning of a period to a new position at the end. 
Crucially, the fluid elements do not necessarily return to the same position — if they did, 
this would be a very poor stirring protocol indeed! 

A periodic stirring protocol in a two-dimensional flow thus induces a homeomorphism if 
from a surface S to itself. A homeomorphism is an invertible continuous map whose inverse 
is also continuous. In our case, describes the mapping of fluid elements after one full 
period of stirring, obtained from solving the Stokes equation, and § is the disc with holes 
(or punctures) in it, corresponding to rods. We treat rods as infinitesimal punctures (see 



Fig. 4(c)), since topologically this makes no difference. As a special case, rods that remain 
fixed are often called baffles. Topologically speaking, moving rods and fixed baffles are 
the same: they are holes in the surface S. However, the homeomorphism (p acts on them 
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differently: the moving rods can be permuted, whilst the fixed baffles remain in place. Here, 
we shall not make a distinction between stirring rods and fixed baffles, and refer to both as 
stirring rods. Hence, a few stirring rods may be fixed by some stirring protocols, such as the 
figure-eight protocol in Fig. [I] which fixes two rods (the islands, or ghost rods). The outer 
boundary of the disc is invariant under ip, corresponding to the no-slip boundary condition. 

Our task is to categorise all possible <p that lead to good mixing. This requires defining 
both what we mean by 'categorise' and 'good mixing.' The categorisation will be done 
up to isotopy, which is a way of defining the topological equivalence of homeomorphisms. 
Two homeomorphisms <p and ip are isotopic if ip can be continuously 'reached' from p 
without moving the rods. If we imagine the two-dimensional fluid as a rubber sheet, this 
means that the two configurations attained by the sheet after application of either tp or ip 
are the same, up to deformation of the sheet. In that case, we write <p ~ ip. Obviously, 
topology is unconcerned by hydrodynamic details and only deals with coarse properties of 
the homeomorphisms. 

To pursue the categorisation, we invoke the Thurston-Nielsen (TN) classification theo- 
rem (HE], which describes the range of possible behaviour of a homeomorphism ip. Specifi- 
cally, the theorem says that ip is isotopic to a homeomorphism ip', where p' is either finite- 
order, reducible, or pseudo-Anosov. These three cases give the isotopy class of ip, and ip' 
is called the TN representative of the isotopy class. Finite-order means that <p>' is periodic 
(that is, p' m = identity for some integer m > 0), and this cannot give good mixing since 
it implies nearby fluid elements will periodically come back near each other. The reducible 
case implies there are regions of fluid that remain invariant under ip', which again is terrible 
for mixing since there are then regions that do not mix with each other. For both finite- 
order and reducible (p f , the actual homeomorphism (p could in practice exhibit much more 
complicated properties than <p', but this cannot be inferred by the rod motion itself. In the 
viscous flows that we have studied, <p and p>' appear to have similar properties. 

The third case, when p>' is pseudo-Anosov (or pA), is both the most interesting mathe- 
matically and most relevant for mixing. In fact, we will define 'good mixing' as <p' having 
the pA property. Mathematically, a pA homeomorphism p>' leaves invariant a transverse 
pair of measured singular foliations, (3^,^) and (3 rs , y u s ), such that <//(3" u , /i u ) = (3" u ,A/i u ) 
and p)'(3 rs ,fi s ) = (3*, A" 1 //)) f° r dilatation A > 1. (The logarithm of the dilatation is the 
topological entropy.) Here /i u and // are the transverse measures for their respective folia- 
tion. There are several terms that need explaining in this definition, and we will illustrate 
what they mean by an example. 

Figure [3] shows the result of a numerical simulation of a two-dimensional viscous (Stokes) 
flow. The container is circular, and the fluid is stirred with four rods, shown aligned hori- 
zontally in the centre. (The velocity field for these simulations was determined using a fast 
spectrally-accurate complex variable method [UJ , and the particle advection computed with 
a high-order Runge-Kutta scheme.) The inset illustrates the motion of the stirring rods: 
they are successively interchanged with their neighbour, in the direction shown. This proto- 
col is written crio^crj 1 , where <7j denotes the clockwise interchange of rod i with rod (i + 1), 
and a™ 1 its anticlockwise counterpart [HI [12J, [19]. Thus, a^ 1 is the inverse operation to crj. 
The subscript i in af 1 refers to the physical position of a rod from left to right, and does 
not label a specific rod. The collection af 1 , i — 1, . . . , n — 1, generates the braid group on n 
strands. We apply the af 1 operations, called braid group generators, in temporal order from 
left to right. Repeated generators are written as powers, as in <7jCXj = of. This sequence of 
generators, called a braid word, defines our stirring protocol, which gives us the homeomor- 
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FIG. 3: Numerical simulation of the four-rod periodic stirring protocol 010203" in a viscous flow, 
corresponding to three rod interchanges for each period (inset). This sequence of interchanges 
forces the flow to be isotopic to a pseudo-Anosov homeomorphism. A material line advected 
for six periods reveals leaves of the unstable foliation 3™. The foliation exhibits four 1-pronged 
singularities around the rods, one interior 3-pronged singularity (in the region indicated by a dot, 
between and just above the first and second rods), and a separatrix attached to the disc's outer 
boundary, visible as an injection cusp along from the top into the mixing region. (The position 
of the singularities and separatrix are approximate.) Compare with Fig. |4j which shows typical 
leaves of the foliation in the neighbourhood of singularities. 



TABLE I: The pseudo-Anosov stirring protocols used in this paper, defined in terms of braid group 
generators. The generator Oj denotes the clockwise interchange of rod i with rod (i + 1), and o^ 1 
its anticlockwise counterpart, where i indicates the physical position of a rod in the sequence from 
left to right. The generators are read from left to right in time. 
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phism ip after we solve the Stokes equations for a viscous fluid. Table [I] defines the protocols 
used in this paper in terms of braid group generators. 

The homeomorphism (p, by the Thurston-Nielsen theorem, is isotopic to the TN rep- 
resentative if', which is pseudo-Anosov in this case. Remarkably, in viscous flows there is 
often very little visual difference between the action of tp and ip', and we find many features 
of ip' reflected directly in Fig. [3j This need not be the case in general: the dynamics of ip 
is at least as complicated as that of <p', in a precise sense [6j |20j I2T] . but can in practice be 
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FIG. 4: For the unstable foliation S^ 1 and stable foliation 3' s : (a) The neighbourhood of a 3-pronged 
singularity; typical leaves of 3^ are shown as solid lines, dashed lines for those of 3~ s . (b) The 
neighbourhood of a boundary singularity. A separatrix of 3^ emanates from the singularity, and 
separatrices for 5P and 3~ s alternate around the boundary, (c) We regard a rod as an infinitesimal 
point (or puncture) , corresponding to a 1-pronged singularity. The leaves of the foliations fold and 
meet at separatrices attached to the rod. 



considerably more so. 

The folded lines in the background of Fig. [3] are a small material closed loop that was 
evolved for six full periods of the stirring protocol. As the material line is evolved for more 
and more periods, it converges to the unstable foliation, tracing out 3 11 and allowing us to 
visualise it as a striated pattern. (The stable foliation 5P is invisible in such experiments, 
so we shall not have much use for it here.) That the foliation 3" u remains invariant under (p' 
means that at each application of <p' the bundle of lines in Fig. [3] is unchanged, except for 
the lines getting denser: each application multiplies the number of lines in a given region 
by A, the dilatation. (A local count of the line density is what the invariant measure /i u gives 
us.) The unstable foliation 3" u is thus the object that captures the essence of stretching and 
folding in a chaotic flow. 

Finally, the 'pseudo' in pseudo-Anosov is tied to the 'singular' in singular foliation. An 
Anosov homeomorphism, such as Arnold's cat map [22], can only exist on surfaces of zero 
Euler characteristic such as the torus. This is because the torus is a surface on which a 
nonvanishing vector field (in this case the foliation) can be smoothly combed. For other 
surfaces, such as our punctured disc, the best one can do is to comb the foliation and 
leave some singularities. Three such singularities are shown in Fig. |4} The singularities are 
characterised by the number of prongs associated with them. The prongs are separatrices, 
emanating from the singular point, around which the foliation branches. For instance, 
Fig. 4(a) shows a 3-pronged singularity: the leaves of the unstable foliation 3" 11 (solid lines) 
branch around the singularity in a triangular pattern around three separatrices. The dashed 
lines show leaves of the stable foliation, 3" s . Figure 4(b) shows a separatrix of 3^ attached to 
the outer boundary of the disc at a singularity. Unstable and stable separatrices alternate 
around the boundary. Figure 4(c) shows a 1-pronged singularity, which occurs around rods, 
as we will see below. 

The stirring device in Fig. [3] shows a total of six singularities in 3^. Most obviously, there 
is a 1-pronged singularity around each of the four rods, because the leaves are folded around 



each rod and meet at a separatrix, as in Fig. 4(c) Next, there is a 3-pronged singularity, in 
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the region marked with a dot in the picture (this dot is not a rod). That the singularity has 
three separatrices is evident if one tries to extend the foliation near the singularity: three 
bundles of leaves will meet at a point, which must then be a singularity. Finally, there is a 
separatrix connected to the disc's outer boundary at a singularity, as in Fig. |4(b)[ though 
this is not as easy to see. From Fig. 4(b) we expect that a boundary singularity will be 
manifested as a cusp in 3^, a consequence of the separatrix emanating from the boundary, 
and leading to injection of material into the mixing region (Fig. [3]). 

In the pA case, almost every aspect of the Thurston-Nielsen classification theorem can 
thus be identified directly in Fig. [3j Rods and the outer boundary possess singularities, and 
often other (interior) singularities arise in the flow itself, such as the 3-pronged singularity 
indicated by a dot in Fig. [3} We will refer to singularities not associated with rods or 
the outer boundary as interior singularities. Our reference protocol of Fig. [3] demonstrates 
that the unstable foliation 3^ and its singularities embody many important features of the 
underlying flow. We shall thus take the unstable foliation as the central focus of our study. 



III. SINGULARITIES OF THE FOLIATION 

An unstable foliation must satisfy three rules if it is to support a stirring protocol corre- 
sponding to a pseudo-Anosov homeomorphisms 



1. Every stirring rod (or ghost rod, such as the islands in the figure-eight protocol) must 



be enclosed in a 1-pronged singularity (Fig. 4(c)). This is a physical requirement: 1- 
pronged singularities are the mathematical consequence of physical stirring, since the 
unstable foliation wraps around the rod. A rod enclosed in a higher-pronged singularity 
makes the rod irrelevant to stirring, since the singularity would exist regardless of the 
rod. Hence, we disregard this possibility. 



2. The outer boundary of the disc contains at least one separatrix of 3" u , as in Fig. 4(b) 
The number of separatrices on the outer boundary corresponds to the number of 
injection cusps into the mixing region. 

3. The smallest number of prongs an interior singularity can have is 3. This is because 
2-pronged singularities are just regular points (they are not 'true' singularities of the 
foliation), and pseudo-Anosovs do not have 1-pronged singularities away from punc- 
tures and boundary components. 

We will use these three rules to limit the number of allowable singularity data of the 
unstable foliation 3 rn . The singularity data of a foliation is the sequence (A sep , N 3 , A 4 , . . .), 
where A sep is the number of separatrices on the outer boundary, and N p is the number of 
interior p-pronged singularities for each p > 3. To enumerate the possible distinct singularity 
data of the unstable foliation for n rods, we use a standard index formula, which relates the 
nature of the singularities to a topological invariant, the Euler characteristic of the disc, 
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FIG. 5: Numerical simulations of (a) the stirring protocol o\a 2 1 V3& 2 1 for 4 rods, with two injection 
cusps; (b) the stirring protocol uxo^ o-&a\ 1 a^a 2 1 for 5 rods, with three injection cusps. 

Xdisc = l. 1 This index formula says that 

n-N scp -J2(p-2)N P = 2 Xdisc = 2. (1) 

p>3 

Formula ([!]) is a well-known extension to foliations of the classical Euler-Poincare-Hopf 
formula for vector fields j25l |26] (see for example page 1352 of (2H p. 1352]). It is important 
to note that the only positive contribution to the left hand side of ([I]) is the term n, the 
number of rods. Hence if there is a large number of rods, there must either be a large number 
of singularities and boundary separatrices, or a small number of singularities some of which 
have many prongs. 

In order for rtlj to be satisfied, the boundary separatrices and interior singularities must 
contribute (2 — n) to the left hand side. For n = 3, the only possibility is to have a single 
separatrix on the boundary. There is therefore unique singularity data for three stirring 
rods in a disc, corresponding to the kidney-shaped region in Fig. [TJ the separatrix on the 
boundary corresponds to a single injection cusp into the mixing region. 

For n = 4, the boundary separatrices and interior singularities must contribute —2 to 
the left hand side of ([I]), which means that there must either be two separatrices on the 
boundary, or one separatrix on the boundary and an interior 3-pronged singularity. These 
two cases correspond to the unstable foliations of the second and third protocols in Table [TJ, 



as depicted in Figs. 5(a) and 3 



1 More generally, the Euler characteristic of a surface of genus g with b boundaries is 2 — 2g — b, where 
the genus is the number of 'handles' attached to a, sphere. In particular, Xspherc — 2, Xdisc — 

2-1 = 1, 

and Xtorus = 2-2-1=0. 
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TABLE II: The allowable singularity data (iV sep , jV 3 , N 4 , . 


. .) for n rods. Each rod has a 


1-pronged 


singularity, and formula ([T| must be satisfied. 


iVgep gives 


the number of separatrices on 


the outer 


boundary, and N p gives the number of p-pronj 


;ed interior 


singularities. 
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6 1 
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6 1 
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6 1 








1 



As we add more rods, there are more possibilities for the singularity data. For example, 



Fig. 5(b) shows a protocol for five rods (n = 5) with three boundary separatrices (A^p = 3), 
corresponding to three injection cusps, and no interior singularities. 

The maximum number of interior singularities occurs when there is a single boundary 
separatrix. The rods and this separatrix then contribute n — 1 to the left hand side of ([!]), 
and so Q can be satisfied if there are n — 3 interior 3-pronged singularities. That is, the 
maximum possible number of interior singularities occurs for N 3 = n — 3. On the other hand, 
in order to have no interior singularities it is necessary to have n — 2 boundary separatrices. 
This corresponds to the maximum possible number of injection cusps into the stirring region. 

Table [TT] provides a complete summary of the allowable singularity data for the first few 
values of n. The number of possible distinct singularity data for a foliation increases sharply 
with the number of rods n, as is evident in Fig. [6] (solid line). A convenient expression for 
the number of singularity data is 



n-3 



# of singularity data = p(k) (2) 



fc=0 



where p(k) is a partition function The partition function counts how many distinct 

ways positive integers can sum to k, with p(0) defined as 1: 

p(k) = # of elements in the set js C Z + : i = fcj. (3) 

The partition function has no simple exact closed form. To find the asymptotic form of ^ 
for large n, we can use the Hardy-Ramanujan asymptotic form for p(k), and replace the 
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10 20 30 

# of punctures, n 



50 



FIG. 6: The number of distinct singularity data increases rapidly with the number of punctures 
(or rods). The solid line is exact and is obtained by summing partition functions as in the 
dashed line is the asymptotic form Q. 



sum by an integral, to get 

# of singularity data 



27^2 (n-3) 



exp (7^2(71 -3)/3) , n > 1. 



(4) 



The dashed line in Fig. [6] shows that the aymptotic form captures the correct order of 
magnitude for large n. 

IV. HYPERBOLIC INJECTION CUSPS 

Injection cusps are not always as plainly visible as in the cases presented thus far. How- 
ever, as we will see in this section, their presence can still be inferred by examining the 
topological properties of the rod motion, including if necessary the motion of ghost rods. 
This section also helps to clarify the type of topological information that can be gleaned 
from the motion of rods [3 [151 128] : 

• If the motion of the rods themselves forms a pseudo-Anosov braid, then the rod motion 
yields interesting topological information even with no knowledge of ghost rods. This 
is the case with all the protocols in the paper thus far (except the figure-eight, which 
involves two obvious ghost rods). The topological information is then "robust," in the 
sense that it does not depend on the specific hydrodynamics of the fluid. 

• If the motion of the rods does not imply a pseudo-Anosov, such as when there is only 
one or two rods [6], and chaotic behaviour is observed regardless, then ghost rods must 
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FIG. 7: (a) Numerical simulation of a single rod tracing out an epitrochoidal path [solid line in 
(b)] stretching a material line for 7 periods. The period-3 orbit discussed in the text is shown 
superimposed, with iterates io, i\, and i% . (b) A Poincare section (stroboscopic map) of some 
representative trajectories shows that the phase space consists of a large chaotic region and several 
regions of regular behaviour. 



be included. This means either looking for regular islands (as in the example in this 
section) or for unstable periodic orbits. However, the presence and location of such 
periodic structures depend on the specific hydrodynamic model (here Stokes flow for 
a viscous fluid). 



The protocol discussed in the present section is of the latter type. Figure 7(a) shows a 
material line advected by a one-rod stirring device, where the rod follows an epitrochoidal 
path. The path of the rod is shown as a solid line in Fig. 7(b)| superimposed on a Poincare 



section. We make two observations about Fig. [7j (i) The Poincare section reveals a mixed 
phase space, consisting of a large chaotic region and several smaller regular regions, including 
a regular region that completely encloses the wall, (ii) The injection cusps into the mixing 
region are not readily apparent, though small cusps are visible. 

The presence of the chaotic region can be understood by examining the motion of the 



physical rod and of the regular islands visible in Fig. 7(b) These regular islands are the ghost 
rods that we use to explain the topological properties of the homeomorphism ip induced by 
the rod motion, as analyzed in Ref. [7]. The braid formed by the rod and the ghost rods 
is a\020zo\<j\a^a2(yz(y\ , which can be shown (with software [29], say) to correspond to a 
pseudo-Anosov isotopy class with a single separatrix on the boundary, as well as two interior 
3-pronged singularities. The singularity data is thus (N sep = 1, jV 3 = 2) for n = 5 rods - 
see Table [TTJ 

This brings us to our second observation: where is the injection cusp associated with the 
boundary separatrix, as predicted by the braid? The injection cusp is there, though it is 
much less evident than those in Figs. [T], [3j and [5j This is because here the separatrix is 
associated with a hyperbolic fixed point, as opposed to parabolic in the previous cases. An 
injection cusp near a parabolic point on the boundary has very 'slow' dynamics near the 
separatrix, meaning that fluid approaches the separatrix very slowly pp. This makes the 
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FIG. 8: Portions of the unstable manifold of each iterate of the period-3 hyperbolic orbit near the 
boundary of the mixing region. The unstable manifold of iterate 12 enters the mixing region and 
corresponds to the location of the injection cusp. 

separatrices clearly visible as unmixed 'tongues' in Figs. [TJ |3j and |5j 

In contrast, the injection into the mixing region is governed here by a period-3 hyper- 
bolic orbit near the boundary of the central mixing region. The orbit — consisting of the 



iterates i , ii, and i 2 — is shown in Fig. 7(a) Notice that the orbit itself does not enter the 



central mixing region, just as in the cases previously considered the parabolic fixed points at 
the wall remain there. However, a portion of the unstable manifold of each iterate is shown 
in Fig. [HJ notice how the unstable manifold of iterate 12 enters the heart of the mixing 
region, following the rod. This unstable manifold is the boundary separatrix predicted by 
the braid. Thus, the way in which fluid enters the mixing region is by coming near 12 and 
then being dragged along its unstable manifold. The iterates io and i\ play no role as far 
as injection into the mixing region is concerned. With hindsight, we can see the unstable 



manifold of i% in the wake of the rod in Fig. 7(a) Because the dynamics in their vicinity is 



exponential rather than algebraic, hyperbolic injection cusps can dramatically speed up the 
rate of mixing [SU] in the central region. The price to pay is an unmixed region around the 
wall of the device. 



V. DISCUSSION 

A stirring device consisting of moving rods undergoing periodic motion induces a mapping 
of the fluid domain to itself. This mapping can be regarded as a homeomorphism of a 
punctured surface to itself, where the punctures mimic the moving rods. Having the rods 
undergo a complex braiding motion guarantees a minimal amount of topological entropy, 
where by complex we mean that the isotopy class associated with the braid is pseudo-Anosov. 
The topological entropy is itself a lower bound on the rate of stretching of material lines, a 
quantity which is important for chaotic mixing. 
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Topological considerations also predict the nature of the injection of unmixed material 
into the central mixing region. The number of boundary separatrices in the pseudo-Anosov 
homeomorphism's unstable foliation determines the number of such injection cusps. The 
number and position of injection cusps is particularly important for open flows, such as 
flows in channels, since there the nature of injection has a profound impact on the shape of 
the downstream mixing pattern (Fig. [2]). 

Topological index formulas allow us to predict the possible types of unstable foliations 
that can occur for a fixed number of rods. We did not provide a way of deriving the 
topological type of the unstable foliation for a given rod stirring protocol. This can be done, 
for instance, by using an implementation [29J of the Bestvina-Handel algorithm [31]. Using 
the enumeration presented here, a mixing device can be designed with a specific number 
of injection cusps into the mixing region, by allowing for enough rods and choosing the 
appropriate stirring protocol. 

More generally, instead of physical rods we can consider periodic orbits associated with a 
stirring protocol. We call such periodic orbits 'ghost rods' when they play a similar role to 
physical rods (that is, material lines fold around them as if they were rods) [7J H2J [15], [28] [32] . 
The topological description presented here applies to the unstable foliation associated with 
periodic orbits. 

In future work, we will consider not just the number of injection cusps, but their relative 
position as well. Indeed, observe that in Fig 5(b) the position of the injection cusps alternates 
sides relative to the line of rods. This is thought to be a general feature of pseudo-Anosov 
stirring protocols, but the proof of this requires careful consideration of whether or not given 
foliations are dynamically allowable, in the sense that they can be realised as the unstable 
foliation of a pseudo-Anosov homeomorphism. 
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